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Abstract 

Kinematical aspects of pion decay tt — > \xv is studied, with neutrino mixing 
^> \ taken into account. An attempt is made to derive the transition probability 

for such a sequence of processes: a ir + produced at (x^^t^) with momentum 

■ p n decays into a /i + and a somewhere in space-time and then the /i + is 
t— I . detected at (x^, t^) with momentum p^ and a u a (a neutrino with flavor a = e, 

fi, ■ ■ • ) is detected at (x u ,t u ) with momentum p v . It is shown that 

, 1. if all the particles involved are treated as plane- waves, that is, if each 

particle is assumed to possess a strictly fixed momentum, the energy- 
momentum conservation would eliminate the neutrino oscillating terms, 
D . leaving each mass-eigenstate to contribute separately to the transition 

^Pi ' probability; 

• !-h ! 2. if one, taking into account that the momenta relevant may not be free 

^ ' from some uncertainty (or dispersion), treats all the particles involved 

■ as wave-packets, the neutrino oscillating terms would appear and would 
be multiplied by two suppression factors, which result from distinction 
in velocity and in energy between the two interfering neutrino mass- 
eigenstates; 

3. as a 2 = a 2 + a 2 + a 2 ,, a n , and a u being uncertainty associated re- 
spectively with p w , and p u , becomes larger (smaller), the feature that 
each of the particles involved propagates along its classical trajectory 
(energies and momenta of the particles involved are conserved during 
the decay) becomes more prominent; and 

4. in the limit of a 2 — > oo, the oscillating terms would again be suppressed 
away. 

An approximate treatment which takes account of the two complementary 
features mentioned above is proposed and similarity and difference between 
our approach and that of Dolgov et al. are discussed. 
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1. Introduction 

Neutrino oscillation [I] is now one of the most exciting topics in particle physics, 
in non-accelerator as well as accelerator high energy physics and in astrophysics [2]. 
There are now plenty of evidence in favor of presence of neutrino oscillations from 
atmospheric neutrino experiments, solar neutrino experiments, reactor neutrino ex- 
periments and accelerator neutrino experiments [3]. 

The standard formulas used to analyze the probability of neutrino oscillation 
between two flavor-eigenstates, e.g., u e and z/ M , read 



where E v and D are the energy and traveling distance of the neutrinos, mi and m.2 
are the masses of two interfering mass-eigenstates, and 9 is the mixing angle. 



is often referred to as the oscillation length. These formulas are usually derived on 
the basis of a plane-wave approach and without paying particular attention to how 
neutrinos are created/detected. 

As is well known, the main physics involved in neutrino oscillation is quan- 
tum mechanics and this aspect of neutrino oscillation has widely been discussed 
in the literature [1] since the pioneering work by Kayser[5J. There are also not a 
few papers which treat neutrino oscillation field-theoretically [6J. In particular, a 
field-theoretical approach of neutrino oscillation with both production process and 
detection process taken into account were recently developed by Giunti et al. in 
[7] and by Asahara et al. in In these works, the neutrino is treated as an 
intermediate particle, while all the external particles are treated as wave-packets. 

We also have, in a series of papers [9] ~ [12], developed wave-packet treatments 
of neutrino oscillation and addressed ourselves, in particular, to such questions as 
"How do neutrinos propagate?", "Equal energy or equal momentum or else?" and 
"How and why does the factor-of-two paradox arise?". In the present paper, after 
extending our previous treatment given in [11] to a three-dimensional case, we shall 
go one step further to study pion decay ir —>■ fiu, with neutrino mixing taken into 
account and with emphasis placed on formal structure of its transition probability. 

ii — > \w decay with neutrino oscillation taken into account has been investigated 
with emphasis placed on such a question as "Do muons oscillate?" by Dolgov et al. 
in [13]. In their treatment, (x ,to), space-time point where and when a pion decay 
occurs, is specified, each of the particles involved is supposed to propagate along its 
classical trajectory, and energy-momentum conservation is imposed by handle They 

1 They, however, at the same time mention that, in quantum field theory, one includes integra- 
tion over space-time point of decay in the definition of the relevant amplitude, and this integration 
leads to the conservation law. 
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first treat all the particles as plane-waves and then introduce momentum distribution 
for the pion alone. 

The present work is more or less stimulated by the interesting work by Dolgov et 
al.|13j and is organized as follows. In Sec. 2 and in Appendix A, we focus on neutrino 
oscillation without paying attention to how the neutrino is created/detected. The 
neutrinos are treated as wave-packets and a number of comments related to those 
questions mentioned above are given. Pion decay with neutrino mixing taken into 
account is investigated in Sec. 3 and Sec. 4. In Sec. 3, we define and derive the transi- 
tion amplitude and probability for a sequence of processes n + — > — > fi + v a with 
all the particles involved treated as wave-packets. The space-time point of creation 
of 7r + and of detection of /i + and v a is specified, while the space-time point of decay 
(x ,t ) is integrated out. The plane-wave limit and the total transition probability 
are also examined. In Sec. 4, we propose an approximate treatment which takes 
account of the two complemetary features: energies and momenta are conserved on 
the one hand and each of the particles propagates along its classical trajectory on 
the other hand. Section 5 compares our approach with that of Dolgov et al. and 
Section 6 gives some concluding remarks. In Appendix B, some algebra relevant to 
Sec. 4 is given. 

2. Three-dimensional wave-packet treatment of 
neutrino oscillation 

Let \u a ) (a = e, /i, • • ■ ) represent neutrino states associated with the electron, muon, 
which are superpositions of the mass-eigenstates having mass m k (k = 
l,2,-.-)H 



t) = ^ ( U ak \u k ). (2.1) 



\Va, 

k 

Suppose that a of momentum p v with uncertainty (or dispersion) a v is created 
at space-time point (x ,t ). Then, its state vector at (x v ,t„) may be written as 

\v^(x v ,t u ;x ,t ;p u ,a u )) = }^ U^kWk)<Pk{x v o, t v o',Pu, (2.2) 

k 

<f> k {x, t; p u , a v ) = (-^) 3 / 2 exp{z(p„f - E k (p u )t) - \o 2 v (x - v k (p u )t) 2 }, 

>7T I 



(2.3) 



where 



X uO % v Xq, t u Q t u tg, 

dE k (p) p 



E k(p) = VP 2 + "I*, v k (p) 



2 We shall assume ^ mi for k ^ I. Also, since we are not interested in CP or T violation, 
we shall take U = (U ak ) to be a real orthogonal matrix. In the two-generation case, U ak may be 
expressed in terms of a single mixing angle 6: U e \ = U^2 = cos 8, U e % = —U^x = sin0. 
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As is well known [H], the wave function of the form given by Eq.(2.3) follows 
readily, if one 

1. starts from a superposition of plane- wave functions!^ 

f°° d 3 p 

(f>k(x,t;p v ,<r v ) = J , 2 , 3/2 f(p;Pv, <r v ) exp{i(px - E k (p)t)}; (2.4) 

2. takes, as the momentum distribution function f(p;Pu,&v), 
f(p;p v ,a u ) = 7 / _^ )3/2 exp{- (P 2 ^' /) }; (2.5) 



3. expands E k (p) around p v as 

E k (p) = E k (p v )+v k (p„)(p-p u ); (2.6) 

4. and performs the p-integration involved in Eq.(2.4) explicitly. 
Note that we have normalized f(p;Pu,cr u ) and (f> k (x,t;p u ,a u ) as 

OO POO 

d 3 p \f{p;p u ,a u )\ 2 = / d 3 x \4> k {x,t;p u ,a u )\ 2 = 1, (2.7) 



and that, to keep this normalization condition valid, we shall define and express the 
plane- wave limit of (p k (x,t;p u ,cr u ) as 

lim (j) k (x,t;p v ,a u ) = — jjt exp{i(p u x - E k (p u )t)}, (2.8) 

that is, letting o v — > in the exponent and, at the same time, on confining the 
neutrino plane-waves within a spatial cube of volume V, letting the normalization 
constant N v = (a v /^f/ 2 -> l/V 1 ' 2 E 

The amplitude for a created at (x ,t ) to be detected as a v a at (^cv, iu) is 
calculated as 

Avn-+v a \ x vOi tuo) = (^a\^fj,(x u , t v \ Xq, t ] p u , O u )) 

= /]Ufj,kUgk < Pk{^vQi'tvG]PvT^v) 
k 

= N u ^ U^ k U ak exp{i(p u x u0 - E k (p u )t u0 ) - ~(?1{x v q - v k (p u )t u0 ) 2 } , 
k 

(2.9) 



3 Neglecting spin degree of freedom, we shall assume that the neutrinos (and the charged leptons 
as well) obey the Klein-Gordon equation. 

4 From the uncertainty relation, one may interpret Kvave-packet = (V^/ a ») 3 as the volume of a 
spatial region within which the neutrino wave-packets are appreciable and express N v — > l/V 1 / 2 

aS Vv^ave — packet y V • 
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and the corresponding propability is calculated as 

Pu^u a (x u0 ,t v0 ) = \A v ^ Ua (x u0 ,t u0 )\ 2 

= ^ UnkUgkUyiUgi (pk(x u0 ,t u0 ]p u ,a u )<pi(x U Q,t U Q-,p u ,cr u ) 



kj 



N l $Z u ^U ak U^U a i exp{-iE [kl] (p v )t v0 } 

k,l 

I I 

e w{-^l(x u0 - v k t u0 ) 2 - -<jI(x„o ~ viUo) 2 } 

k,l 

exp{-al(xvo - v { ki)(Pu)tuo) 2 }exp{-^al(v M (p v )) 2 tl }, (2.10) 



where 



E {k i){p u ) = ^{E k (p v )+Ei(p v )), E [kl] (p u ) = E k (p v ) - Ei(p v ), 

V(M)(Pu) = ^(v k (Pu) +vi(pu)), v [kl] (p u ) = v k (p v ) - vi(p v ). 
With Eq.(2.7), we have 

[ d 3 x u P u ^ Ua (x u0 ,t u0 ) = 1. (2.11) 
A couple of remarks are in order as regards Eq.(2.10). 

1. Each term in Eq.(2.10) contains explicitly a factor which implies that the 
space-time point (x v ,t v ) where and when the neutrinos are detected nearly 
satisfies 

x v = xo + V(ki)(p v )(U-t ), (2.12) 

which is precisely what could be referred to as classical trajectory of the neu- 
trinos. We may therefore, assuming that (x u0 — V( k i)(p v )t v0 ) 2 <C cr~ 2 is satisfied 
for any k and /, express Eq.(2.10) approximately as 

- N l Yl u ^U ak U^U a i cos(E [kl] (p v )t v0 ) exp{~al(v [kl] (p u )) 2 tl } 
k,i 

AT 2\^TT TT TT TT (n D \ ( 1 2 (%i](P")) 2 n 2l 

= N; > U^UakU^Uai cos(2tt- ) cxp{--a v 1 ' U2 J }, 

(2.13) 

where D = \x v q\ = \v(ki)(Pv)\t v o is the traveling distance of the neutrinos, and 

27r|ff (fe0 (p l/ )| 4nE (kl) (p v )\v {kl) (p u )\ 
\E[ki]{Pu)\ 



\m% - m. 
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2. In each of interference terms, the oscillating factor cos(E^i](p u )t u0 ) is multi- 
plied by a factor, denoted by which implies that this term is significant 
only when the condition 

t„-t <-J—: or D<*p^, (2.15) 

CTv\V[kl]{Pv)\ CT u \V[kl]{Pu)\ 

is satisfied. This effect, often referred to as coherent condition for neutrino 
oscillation [15], manifests itself explicitly in our wave-packet treatment but not 
in conventional plane- wave treatments. Note also that Vk{Pu) 7^ vi{Pu) for 
k I only results in appearence of this effect, but by no means affects the 
oscillating factor and hence the oscillation period 2%/ E[ki](p u ) or oscillation 
length l k ^iE 

3. In conventional plane-wave approaches, some prefer to work with an "equal- 
momentum prescription" and others an "equal-energy prescription" [5] . [16] . 
We like to point out that the wave-packet treatment _given here reduces smoothly, 
in the plane- wave limit, to the former prescription.] and that, although these 
two prescriptions give approximately the same oscillation formulas for rela- 
tivistic neutrinos, they have to be distinguished from each other in general 
and conceptually [TO] . 



3. Wave-packet treatment of n — > [w Decay 

3.1. 7r — > \iv [L decay followed by — > v a transition 

We now attempt to formulate such a sequence of processes: a 7r + produced at (x n , t n ) 
with momentum p n decays into a fi + and a somewhere in space-time and then 
the fi + is detected at with momentum p^ and a u a is detected at (x u ,t u ) 

with momentum p v . Taking into account that the momenta relevant may not be 
free from some uncertainty or dispersion, we shall treat all the particles involved as 
wave-packets and examine the plane-wave limit later on. 

If the decay space-time point is denoted by (xo, t ), we may use the amplitude we 
have introduced, Eq.(2.9), to describe a created at the decay point to be detected 
as a v a at (x u ,t u ). Similarly, the amplitude for a ji + created at the decay point to 

5 If, noting that the magnitude of the amplitude A v ^^(x^o,^), Eq.(2.9), has a peak at 
= Vk{pv)t V Qi one substitutes this into its phase factor to obtain exp{— i(myEk(pv))t v o}, the 

phase factor of the probability (2.10) would become exp{— i{m\/ Ek{p v ) — mf j Ei(p v ))t u o} , which 
will results in an oscillation period nearly twice as large as the standard one for relativistic neutri- 
nos. This is where the so-called " factor-of-two paradox" [16] arises and is obviously an inadequate 
prescription. 

6 We shall give a more general three-dimensional wave-packet treatment of neutrino oscillation 
in Appendix A. Note that, for relativistic neutrinos and in the two-generation case, Equations 
(2.13) and (2.14) coincide (apart from the normalization factor) with Eqs.(l.l) and (1.2), if one 
let a v =>■ and identifies Er^n^u) with E u . 
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propagate to its detection point (x^t^) and the amplitude for a tc + produced at 
(xtt^t,-) to propagate to its decay point may be described respectively b}0 

A^x^o-, tfio) = 0/X x mo> ^o>2V> °/J 

d 3 p 



^ 27r y /2 f(F,P», <7») exp{i(px M o - E^p)t^)} 
= exp{i(p M f M o - Eptpo) - ^jfeo - v^t^) 2 }, 

A^XQtt-i to-ir) = 0vr(a;O7r, to-*] Pit-, C"tt) = ^(^irO) ^ttO! Ptt, 0"tt) 

, 2 s 3/2 f(p;Pir, exp{-i(pf^ - K(p)^ )} 

= iV 7r exp{-i(p 7r ^7rO - KUo) - ^l(x n0 - ^tto) 2 }, (3.1] 
where and a n are uncertainties associated respectively with p^ and and 

^A = ("^) 3/2 , ^AO — 2?A X , ^A0 = t\~ t , 



'It 

E\ = Jjfx + ml, v x = tt, A = 7r. 

The transition amplitude for 7r + — > /U + z/ M — > yU + z/ a with the decay space-time 
point (xo,to) specified is then given by 

A^— » / Lti/ p — >fj,u a \ X 0i ^o) = A v ^ Vol {x v 0i tuo)A^(x^Q, t^An^Xoir, to w ) 

= y^^UfMk UgkAu (xp ,t ), (3.2) 
k 

where 

Ak(x ,t ) = 0fc(^O,^O;^,^)0A t (^O,^O;P At ,^)C(^7rO^7rO;P7r,CT7r) 

= N u exp{i(p u x u0 - E k t u o) - -&l{xvo - VkUo) 2 }, 

N n eocp{-i(^ r x w o - E^o) - ^(^o - ^Ao) 2 }- (3.3) 

The transition amplitude for the whole process we have specified in the begining 
of this subsection may be obtained by integrating Eq.(3.2) with respect to x and 



>/m^— >/j,v a — J d xo J dto A ir —fi U , ll — > .f iVa {xo,to) 



7 From now on, momentum dependence of energies and velocities will not be indicated explicitly, 
if not particularly necessary. 

8 It is understood that the ^-integration extends over a infinitely large time interval T. 
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and the corresponding transition probability per unit time interval is given by 

= 'y^^UfikUakUfdUal Gkh (3-5) 

k,l 

where 



A k = J d x J dt A k (xo,t ), (3.6) 
G k i = 7pA k A* = ^ J C dx J dt J dx J dt' A k (x o ,t )A*(tf o ,t' o ). 



(3.7) 



3.2. Transition amplitude and probability 

To calculate A k , Eq.(3.6), we write Eq.(3.3) as 

1 



A k (x ,t ) = Nexp{i(6 k - Apx + AE k t ) - ^a 2 y k (x ,t )}, (3.8) 



where 



7TV7T 



Ap = Vu+Vfi- Pn, AEfc = E k + E^- E v , 
k Pv%v E k t u -\- Pf^Xfj, E^tfi Ptt-Ett E w t w , 

yk{x , t ) = ^-{al(x u0 - VkUof + cr 2 (x^ - v^t^) 2 + o\{x^ - v^Uo) 2 }. 
a 



(3.9) 



Introducing^] 



X k X V V k t Ul X n Vutif, X-ft X w V^t 



(u) k = (alu k + <rJu M + alu n )/a 2 , 
{uw) k = {a 2 v u k w k + alu fl w tl + alu n w n )/a 2 
(u,w) k = (uw) k - (u) k (w) k , u,w = v,X, (3.11) 



9 (u, w) k may be expressed also as 

1 



(u,w)k = +^^5[k]»H +^^V1'"m}' ( 3 - 10 ) 

where 

U[ K \] = U K - U\, K, A = k, /i, 7T. 
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we further express yk(xo,t ) as 

Vk{xo,to) = (x - (v) k t - (X) k ) 2 + a k (t - t k ) 2 + c fc - a k t 2 k , (3.12) 
and i(0 k - Apx + AE k t ) - \(y 2 y k {x ,t Q ) as 

i(9 k - Apf + AE k t ) - icr 2 y fc (^o,^o) 



i(0 fc - Ap(X) fc + AE k t 



^cr 2 (f - C^Vo - (^)fc + i-zAp) 2 - -a 2 a k (t -t k - i— — AE k ) 2 , 
2 a z 2 o z a k 

(3.13) 



where 



a k = (v,v) k , b k = (v,X) k , c k = (X,X) k , 

t k = -h/a k , AE k = AE k -(v) k Ap. (3.14) 

Substituting Eqs.(3.8) and (3.13) into Eq.(3.6) and performing the x - and t Q - 
integrations over the whole space-time, we find 

lix 9tt 

A k = N(—f 2 (^)^ 2 exp{t(e k -Ap(X) k + AE k t k )} 

<7 Z <J CLk 

exp{-^(Ap) 2 - 1 ±-(AE k ) 2 - l -a 2 {c k - a k t 2 )}. (3.15) 
Substituting Eq.(3.15) into Eq.(3.7), we obtai 
G ki = — A k A* 

1 ?7T 9TT — 

= - N {—f 2 (— ) 1/2 exp{i(6 k - Ap(X) k + AE k t k )} 
1 a z <J z a k 

-P{-^(APT - - \^c t - a k tl)) 

2ir 2ir -> 

N Hf 2 (~^) 1 ' 2 exp{-i(0, - Ap(X) l + Am)} 
a z & z ai 

= N H exp{-a 2 Z M - \h h - iQ k i}, (3-16) 

a z 



10 One may substitute Eqs.(3.8) and (3.13) directly into Eq.(3.7) and perform the Xq-, to~, x'q- 
and ^-integrations, by changing the integration variables from xq and x' to x_ = Xq — Xq and 
x + = (xq + x' a )/2, and from to and t' to i_ = to — t' a and t + = (a k to + a*it'o)/(a k + a;), to obtain 
the same result. 
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where 




(3.17) 



(3.20) 



(3.19) 



(3.18) 



2. 



couple of remarks as regards "trajectories of particles" are in order. 

From Eq.(3.8), one sees that \Ak(xo,to)\ has a peak at yk(%o,to) = 0, which, 
from Eqs.(3.9) and (3.12), implies on the one hand that 

Xvo - v k Uo = 0, x m0 - v^t^o = 0, x n0 - VirUo = 0, (3.21) 

and on the other hand that 

x - (v)kto - (X)k = 0, t - t k = 0, c k - a k t 2 k = 0. (3.22) 

The three equations in Eq.(3.22) may be derived somehow directly as some 



weighted averages of the three equations in Eq.(3.21)0 

From Eqs.(3.15) and (3.16), one sees that \Ak\, as well as Gkk, has a peak at 
Cfc — o-k^k — 0- Thus, even after the integrations with respect to x and t have 
been carried out, the third equation in Eq.(3.22) remains meaningful and may 
be regarded as something which reflects that each of the particles involved in 
the decay propagates along its classical trajectory, Eq.(3.21). 

For k /, Ck — akt\ = and c\ — aitf = are not compatible with each other 
and Zki = (cfc — akt\ + q — a^)/2 = never holds. Thus, for the interference 
terms, Gki with k ^ I, it is actually not quite clear as to if and how trajectories 
may be defined for the particles involved in the decay. We shall come back to 
this question in the next section. 

lote that the three equations in Eq.(3.21) may be rewritten as 

xo-Vkto = X k , xq-v^q = Xp, x - vUo = X-x, (3.23) 

which it follows that 

v\k^]ta = -X[kn], v^]to = - X^], w[7rfe]io = — X^k]- (3.24) 

irst equation in Eq.(3.22) follows from Eq.(3.23), while the second and third equations in 
.22) follow from Eq.(3.24). 
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3.3. Plane-wave limit 



Let us examine the plane-wave limit: a — ► or, eqivalently, ov, cr M , cr^ — > 0. By 
going back to Eqs.(3.3) and (3.6) and trivially performing integrations with respect 
to xq and to, one finds 

lim A k = -4- (2vr) 4 exp{z0 fc } <5 3 (Ap) 5(AE k ), 



from which it follows that] 12 ! 

lim G w = I -L (2tt) 8 exp{-^ N ^} [<5 3 (Ap)] 2 5(A£ fc ) «5(AS,) 

= 5 W ^ (2vr) 4 <5 3 (A$ (J(AS fc ). (3.25) 
Substituting Eq.(3.25) into Eq.(3.5), one obtains 

P^»^» a = ^ (2tt) 4 £ U 2 k U 2 ak 5\p u + p, - p.) 5(E k + E,- E n ), 



or, in the two-generation case, 

P*-*^-*^ = (2vr) 4 5 3 {p u +p^- p n ) 



(3.26) 



{sin 4 9 5(E 1 + E^ — + cos 4 6 5(E 2 + E^ — E w )}, 

P^u^ Ve = (2tt) 4 5 3 (p; + p^ - p n ) 

sin 2 6 cos 2 {5(^1 + E„ - E n ) + 5(E 2 + E„- E n )}. 

It is seen that energy-momentum conservation prevents different mass-eigenstates 
to interfere with one another and, as a result, each of mass-eigenstates appears to 
contribute separately to the transition probabilities!* 3 ! 

It is also interesting to examine the cases in which one of a u , and a n is kept 
finite. From Eq.(3.15), with a n kept finite, one finda^ 4 ! 

lim A k = i N„ (^) 3 / 2 (2n)6(AE k - v n Ap) 

exp{i(p u x vn - E k t un +p^x^ n - Efjt^) - ^(Ap) 2 }, 



12 To derive Eq.(3.25), it is taken into account that, for k ^ I, AE k = and AEi = are not 
compatible with each other and use is made of the farmiliar technique to handle with the square 
of a S-i unction 1 17]. 

13 An intuitive argument on connection between neutrino oscillation and energy-momentum 
conservation has been given by Kayser[5] and by Lipkin[T5]. 

14 x K .\ = x K — x\, < k a = t K — t\, k,X= V) jU, 7r. Note that, with 0^,07, — > 0, one has a — > a„, 
AE k — > AE k — v-nAp, (X) k — ► X„ and afc, 6fc, c k — » 0, while t k remains finite. 
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from which one derives 

lim G kl = i 4 N * (^) 3 (^)^( AE k ~ v,Ap)5{AE l - v n Ap) 

eacp{-iEr fc nt w - -\(Ap) 2 } 

= 5 kl i (2vr) 4 (^) 3 / 2 5 ( - v n Ap) exp{-4(Ap) 2 }. 

(3.27) 

It is interesting to see that there appears in the amplitude a single ^-function which 
implies that energies are conserved in the rest frame of the pion and this ^-function 
eliminates interference terms from the transition probability. In contrast, with a v 
kept finite, one has 

lim A k = i N v & 3 / 2 (27i)5(AE k - v k Ap) 



a 



exp{z(p At f (UV - Ept^ -p n x wu + E n U v ) - — ^(Ap) 2 }, 
lim G kl = 4 (2vr) 4 (^) 3/2 ^ «y(AS fc - v k Ap*)5{AE l - v t Ap) 

exp{-^(Ap) 2 }. (3.28) 

Here, remarkably, interference terms remain nonvanishing, but appear not to contain 
an oscillating factor! 1 ^! 

3.4. Total transition probability 

From Eqs.(2.4) and (2.8), noting that the function f(p;p v , a v ), defined by Eq.(2.5), 
satisfies 

d?Pv\f(p;pv,Vv)\ 2 



one may derive such a relation aa 16 l 

d 3 p u / d 3 x u faixvOftriip^aJtftfriJriiPwa,,) 



dp v \ d x v lim .\^k(.x^,t^p v ,a v )<j>i(^,t M ;p v ,(T v )]. (3.29) 



15 AE k — v k Ap — and AE[ — ViAp — are not incompatible with each other and the product 
of the two (5-functions, 6(AE k — VkAp)8(AEi — viAp), may be rewritten as S(AEk — v k Ap)5((Ek — 
Ei)(E k + Ei + E,j, - E n )/Ek). Note also that, if one lets => in Eq.(3.27) or lets a v in 
Eq.(3.28), these equations reduce to Eq.(3.25). 

x '\0 = x * — x 0> ^AO = ^ A — ^0! K = ^ M) 7r - 
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On substituting Eq.(3.3), we integrate Eq.(3.7) with respect to x v and p u and consult 
Eq.(3.29), to obtain 



G 



ki 



J d 3 p u J d 3 x v 

= ^ J d 3 x J dt J d 3 x' J dt' Q ^(fno,^;^,^)^^,*^;^,^) 

M (^o, tuoipp, 0-^)0* (x 1 ^, t^ ; pp, (7 M ) 
d 3 p u d 3 x v lim [^(^0,^0;^/,^)^* (^ ,^ jp,,,^)] 

= / d 3 ^ / d 3 ^ lim (3.30) 
Similarly, we obtain 

J d%J d 3 Xlx G kl = J d 3 p tl J^d 3 x ll Hm o G fcJ , (3.31) 

y d%J d 3 x^ J d 3 p v J d 3 x u G k i = J d 3 p^ J d 3 XfM J d 3 p v J d 3 x v ^hm^G^. 

(3.32) 

Equation (3.32) implies that, as far as the integrated transition probability is con- 
cerned, there is no difference between treating decay products as plane-waves or as 
wave-packets. 

Integrating Eq.(3.5) with respect to x u , p u , x^ andp M , summing it up with respect 
to a and consulting Eqs.(3.32) and (3.27), we obtain the total transition probability 
or total decay rate as 

PjT^flU^ = ^ ^ j d P^ [ d X^ [ d p u [ d X v Pw^flU^^^tVa 

a. 

= ^U 2 k / d 3 Pll / d 3 x^ / d 3 p v / d 3 x v lim G kk 



M^A^E U% [d 3 pjd 3 p u 

™$ v J J 



1 



5(AE k - v n Ap) exp{ -(Ap) 2 }, (3.33) 



which reduces, in the limit of a n — > 0, to 

= (27T) 4 U lk j ^Vv 

$ 3 {Pv + P»~ Pn) S{E k + E^ — E n ). (3.34) 

Equation (3.34) as well as Equation (3.26) are expected results, which serves as a 
check of consistency of our treatment as a whole. 
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4. An approximate treatment 



In this section, we like to propose an approximate prescription to handle with the 
trajectory factor exp{— a 2 Zki}, the energy-momentum factor exp{— H^/a 2 } and the 
phase factor exp{— iOfczjJEjj 

4.1. Trajectory factor 

As noted before, the factor exp{— a 2 Zkk\ = exp{— a 2 [ck — a^)}, contained in Gkk, 
may be interpreted as something which reflects that each of the particles involved 
propagates nearly along its classical trajectory, Eq.(3.21), while, for k ^ I, Ck~akt\ = 
and q — ait 2 = are not compatible with each other and Zki = {ck — a^t 2 . + c\ — 
ait 2 ) /2 = never holds. We have encountered a similar situation in Sec. 2. There, 
writing 



1, , . 

v k: i = v (ki) ± -V[ M] , (4.1; 



and, accordingly, 
1 



((x v0 - v k t u0 ) 2 + (x„ - ViUo) 2 ) = (x u0 - v { ki)t u0 ) 2 + ^(v [H ]) 2 t 2 u0 , 



we have reasonably regarded Eq.(2.12) as a classical trajectory of the interfering 
neutrinos. Here, on writing Zki as 

Zkl = C{kl) — a (kl)tfkl) + (kl, (4.2) 

and, accordingly, exp{— a 2 Zki\ as 

ex.p{-a 2 Z k i} = z H C k ] \ (4.3) 



where 



we shall regard 



z k i = exp{-a 2 (c ik i) ~ a { ki)t 2 {k i))}, 

$ = exp{-a 2 C fc z}, (4.4) 
C(ki) - a {k i)t\ki) = (4.5) 



as something which reflecs that each of the particles involved, including the inter- 
fering neutrinos, propagates along its classical trajectory and hence z k i as a factor 



17 We shall introduce quantities with suffix (kl) and quantities with a prime or a double prime; 
the former are, by definition, quantities depending on v^i) = (wfc+uz)/2 and/or = (Ek+Ei)/2, 
just as the quantities with suffix k depend on Vk and/or E^, while the latter are defined explicitly 
in Appendix B and vanish if k ^ I or in the limit of <j v — > 0. Also, many of equations given in this 
section will be derived explicitly in Appendix B. 
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representing that each of the particles involved propagates nearly along its classical 
trajectory. 

(ki, defined by Eq.(4.2), is given explicitly by 

Cki = -{c k - a k tl + ci - aitf) - (c( k i) - a( k i)t 2 (kl) ) 

1 

a (kl) ((a m + a'^) 2 - (a' kl ) 2 ) 

x[ (a(fci) + a'Diafkifki ~ ^(ki)b( k i)b kl + h \ki) a ki + a (fco( a K - (Kif)} 
~ ( a (ki)b' k i ~ b {k i)a' H ) 2 - a (k i)(c' kl (a' H ) 2 - 1b" kl d kl b' kl + d kl (b' kl ) 2 ) ]. (4.6) 

When Equation (4.5) holds, which implies 

X[ kX] = -V[ K \]t , re, A = (hi), li, 7T, (4.7) 
t being a constant independent of k and A, ( k i reduces to ( k i(t u — t ) 2 , where 

j ( a (ki) + a ki)( a (ki) a ki ~ \( a 'ki) 2 ) ox 

U ' " ( a (kl) + a kl) 2 - ( a 'kl) 2 ' 1 ' ' 

and the trajectory factor exp{— a 2 Z k {\ may be approximated as 

exp{-a 2 Z kl } w zhSm (4.9) 

where 

t£ = exp{-a 2 C kl (t u -t ) 2 }. (4.10) 

Necessary conditions for interference terms to be appreciable follow from Eqs.(4.3), 
(4.4), (4.9) and (4.10): 

(a) C( H ) - a(fcz)%) < I/a 2 ; and 

(b) C« or C«(<,-<o) 2 < V^ 2 . 

The former condition applies to the diagonal terms G kk too and implies that each of 
the particles involved should propagate along its classical trajectory at least approxi- 
mately, or that those mass-eigenstates which do not satisfy this condition contribute 
little to the transition probability, while the latter is a condition which corresponds 
to Eq.(2.15), known as coherent condition for neutrino oscillation [15], and implies 
that those two mass-eigenstates, if their masses do not satisfy this condition, inter- 
fere little with each other. 

From Eq.(4.7), it follows that 

X v - V(H){tu ~ t Q ) = x^-v^t^-to) = x n + v n (t -t n ) = x . (4.11) 

Since these equations appear to be exactly similar to Eq.(3.21) which involves the 
decay space-time point (xo,to), we shall refer to (xo,to) as pseudo decay space-time 
point. Note that the pseudo decay time to, and hence the pseudo traveling time of 
the neutrinos t„ — to as well, may in principle be deduced from knowledge of [x v , t u ) 
and VQsi) and of (x^,t^) and (or of (x w ,t n ) and v w ). 
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4.2. Energy-momentum factor 

In the limit of a — > 0, exp{— H kk /a 2 } gives rise to 6(E), + E^ — E_ B -)6 3 (p u + p^ — 
p w ), while exp{— H k i/a 2 } with k ^ I vanishes (see also Sec.3.3)o In order for 
interference terms to remain nonvanishing and appreciable, the two constituting 
factors of exp{-H kl /a 2 }, exp{-(Ap) 2 /2a 2 -(AE k ) 2 /2a 2 a k } and exp{-(Ap) 2 /2a 2 - 
(AE)) 2 /2a 2 ai}, have to overlap appreciably with each other. 

To look into this condition more deeply, we rewrite (AE k ) 2 /2a k + (AE{) 2 /2ai as 

-L(AE k ) 2 + ^-(AE t ) 2 = J-(AE m ) 2 + Vkl , (4.12) 
la k 2ai a(fcj) 

and, accordingly, exp{— H k i/a 2 } as 



where 



exp{--^} = hu$\ (4.13) 



h kl = exp{-±((Ap) 2 + —(AE {kl) ) 2 )}, 
<? 2 a(ki) 

eg 3 = exp{-l^}. (4.14) 
a z 

The factor h k i has a peak at 

Ap = Pv + Pn — P-k = 0, 
AE m = E (kl) + E^-E^ = 0, (4.15) 

which may reasonably be regarded as representing the energy-momentum conserva- 
tion in the presence of neutrino mixing. 

rj k i, defined by Eq.(4.12), is given explicitly by 

m = ^-(AE k ) 2 + ^-(AEi) 2 - —(AE (kl) ) 2 
la k Zai o^i) 



4a (fcI) ((a (JH) + < 2 ) 2 - (a' kl ) 2 ) 

x[ (a(fci)(o(fci) + a ki)( AE 'ki) 2 ~ 4a( H )a' kl AE {kl) AE kl 
+ 4((a' kl ) 2 - 4(a (w) + 4))(A£ (fc0 ) 2 ]. (4.16) 



18 More precisely, one has 

lim exp{-^H kl } = 6 kl (ira 2 ) 3 / 2 {^o 2 ^) 1 ' 2 S 3 (Ap) S(AE k ), 



cr^O 

and, accordingly. 



lim G w = S kl I l^) 1 / 2 (2^) 4 <S 3 (Ap) <5(A£ fc ). 
In order for this expression to coincide with Eq.(3.25), one needs to identify T with {TT/a 2 a k ) 



1/2 
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When Equation (4.15) holds, r) k i reduces to 



O(fcj) +a kl 2 

and exp{— Hki/cr 2 } may be approximated as 

ex p{-i_iJ H } » ^g), (4.18) 



where 



5? = exp{-l^}. (4.19) 



Necessary conditions for interference terms to be appreciable now follow from Eqs.(4.13), 
(4.14), (4.18) and (4.19): 

(c) (Ap) 2 + (AE {kl) ) 2 /a {kl) < a 2 ; and 

(d) r] k i or rj k i < a 2 . 

The former condition applies to the diagonal terms G kk too and implies that energies 
and momenta involved should conserved at least approximately, or that those mass- 
eigenstates which do not satisfy this condition contribute little to the transition 
probability, while the latter condition implies that those two mass-eigenstates, if 
their masses do not satisfy this condition, interfere little with each other. It is to be 

(2) "(2) 

noted also that Cm or Cm seems to have something to do with the suppression factor, 
given by exp{ — {E^) 2 /2a 2 {yl + vf)}, which appears in a more general wave-packet 
treatment of neutrino oscillation (see Appendix A). 

4.3. Phase factor 

6fcz, Eq.(3.20), may be rewritten as 

e« = AE'uiU-^tk + U))- AE [kl) {t k -U) 

- \w' S+ , ( a (ki) + OkiHhki} + Ki) - a ki b 'ki -i 

(fl(fci) + a kl) ~ KJ 

, 9A P (Ki( a (ki) + a 'Li) ~ a'ki(hki) + b' u ) 

+2AEm — ^JW^W — • (420) 

When Equation (4.7) holds, Q k i reduces to Q k i(t u — to), where 

n a g,/ Q(fcQ(«(fcQ + a'L) ~ (a'ki) 2 / 2 A p (<*(*«) ~ a 'Li) 

U kl = AE kl — — — — AE {H) 



(«(«) + O 2 - Ki) 2 ( ' («(«) + O 2 - Ki) 2 ' 

(4.21) 

If Equations (4.15) are further applied, Q k i reduces to Qki{tv — to), where 

f, ajkiMk^ + a'D-Ki) 2 / 2 ^ (A00 , 

° H = (a m + <i) 2 -K) 2 E[klV (422) 



17 



4.4. Summary 

We have rewritten Eq.(3.16), 

G k i = N M exp{-a 2 Z M - -^H M -i9 M }, 

as 

Gki = N u z u h kl $ $ exp{-t& kl }, (4.23) 

where z k i, h k i, and are given by Eqs.(4.4) and (4.14). Substituting Eqs.(3.16) 
and (4.23) into Eq.(3.5), one obtains 

P-K^ixv^iiva = U^kUptkUixiUgi N ki exp{—<7 2 Z k i -Hki — iQki} 

k,i ° 

= Yl U ^U ak U,iU al N kl z M h kl $ $ cos(e w ), (4.24) 

k,l 

We have then shown that Equation (4.24) may be approximated as 

~ '^^U^kUakU^Uai N k i z k i h k i cos(Q kl(t u — t )), 

k,l 

(4.25) 

where , , and Q k i are given respectively by Eqs.(4.10), (4.19) and (4.22). 

We have furthermore discussed implications of each factor in Eqs.(4.24) and 
(4.25) and derived necessary conditions for oscillating terms to be appreciable. It 
is to be emphasized here that the two features implied respectively by z k \ and h k i, 
that is, each of the particles involved propagates nearly along its classical trajectory 
on the one hand and energies and momenta of the particles involved are nearly 
conserved on the other hand, are complementary to each other in the sense that, 
as a becomes larger (smaller), the former (latter) becomes more prominent, and 
that the two suppression factors ^ or £ kl and ^ or ^ are also complementary 
to each other in the sense that, as a 2 becomes larger (smaller), the former (latter) 
becomes more effective. 

If C(fcj) — a(ki)t 2 k i) <^ l/o" 2 is assumed to hold for any k and /, Equations (4.24) 
and (4.25) may be approximated as 

^ U VkUakU^lU a l N kl h k l £ k ) f£? C0S(6 H (^ - t )) 

k,l 

» U ^U ak U^U al N kl h kl t k i i { k f cos(t fc/ (t, - *o)), (4-26) 
k,i 

where Q k i is given by (4.21). Equations (4.25) and (4.26) are to be contrasted with 
Eqs.(2.10) and (2.13) and also with Eq.(3.26). 
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5. Comparison with the approach by Dolgov et 
al. 



In the wave-packet treatment of n — > fiu decay given in [13J , assuming that a pion 
has been created with some momentum distribution f(p n ), the authors define a 
transition amplitude and the corresponding transition probability by 

A ^nu^u a (^t ) = ^U^ k U ak / d 3 p w f(p w ) exp{iip%(p w )}, 



and 



'y~^ J UfikUakUfdUal G kl (x ,t ), (5.1) 



k.l 



where 



GZ@o, to) = J dPp* J d s p'J(p n )m) eMK^kiPn) ~ </>?(&))}, 



with 



Vk (Ptt) = PukXuQ - E k (p vk )t u Q + p^kX^Q - E^p^t/jjo - PnXnO + E w (p n ) t ^0 , 



Here, the energy and momentum of the neutrino mass-eigenstate and of the muon 
are all regarded as depending, through the energy-momentum conservation law, on 
the mass of the neutrino and hence carrying the suffix k or /. Assuming that the 
dispersion of f(p n ) is small and that each of the particles involved should propagate 
along its classical trajectory (cf. Eq.(3.21)), they claim that Equation (5.1) reduces 
to 

2 2 

P^u^S^o) = ^U^UckU^Uca cos( m * ~ m ' U ). (5.2) 



k.l 



v 



Although we may also define an (x , £ )-dependent transition probability from 
Eqs.(3.2) and (3.3), we like to point out here that, in spite that the approach we 
have developed in Sec. 3 and Sec. 4 is characterized among other things by the space- 
time point of decay, (x , t ), being integrated out, the transition probability we have 
derived, Eq.(4.26), seems to describe a situation very much close to what Dolgov 
et al. claim to describe with Eq.(5.2): Equation (4.26) is derived with each of the 
particles involved assumed to propagate along its classical trajectory, and, as a result, 
becomes depending on the pseudo traveling time t u — to, which just corresponds to 
the traveling time t„ — t in the (xo, to)-dependent probability (5.2) of Dolgov et 
al.. It is to be noted at the same time, however, that Equation (4.26) appears to 
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be distinct from Eq.(5.2) in the following three points: (a) each term in Eq.(4.26) 
contains a factor which implies that the energy-momentum conservation holds only 
approximately; (b) each of the oscillating factors in Eq.(4.26) is multiplied by the two 

suppression factors £^ and which are absent in Eq.(5.2); and (c) the oscillation 
period in Eq.(4.26) deviates in general from the oscillation period 4irE u /\ml — mf\ 
in Eq.(5.2). 

In addition to the case in which both neutrinos and muons are detected (Case 

A) , Dolgov et al. [I3j consider also the case in which only muons are detected (Case 

B) and the case in which only neutrinos are detected (Case C). In our approach, 
the transition probability corresponding to Case B may be obtained by integrating 
Eq.(3.5) with respect to x v and p v and summing it up with respect to a: 



— ^ ^ / d Pu I d x v P-n- 
= ^ UvkUakU^Uai J d 3 p v J d 3 x u G k l 

a,k,l 

'% I d z p u I d 3 x u G kk . (5.3) 



As seen, interference terms drop out, implying that muons do not oscillate. Although 
our approach gives, not only for Case A but also for Case B and Case C, results more 
or less different from what Dolgov et al. claim, we share with them the conclusion 
that muons do not oscillate!^! 



6. Conclusions 

Our wave-packet approach to pion decay is characterized by treating all particles in- 
volved as wave-packets and by integrating out the space-time point of decay (x , t Q ). 
We have seen that, in such an approach, (1) energy-momentum conservation appears 
to hold only approximately, (2) exact energy-momentum conservation, which holds 
in the plane- wave limit, would eliminate neutrino osillating terms from the transition 
probability, and (3) treating only one of the particles involved as a wave-packet is 
insufficient to allow for neutrino oscillating terms to appear. We have furthermore 
developed an approximate treatment, which allows one to define a pseudo decay 
space-time point (x ,t ) and to express the transition probability, originally given 
by Eq.(4.24), approximately as Eqs.(4.25) and (4.26). 

To conclude, we like to mention that phenomenological implications of the out- 
comes of the present study need to be examined separately and that our approach 
may readily be applied to various decay as well as production processes. 



19 Note however that there have been some debates as regards whether or not muons can 
oscillate [T§]. 
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Appendix 

Appendix A. A more general three-dimensional wave-packet 
treatment of neutrino oscillation 

More generally, one may conceive a situation in which the peak momentum of the 
neutrino mass-eigenstates, p v in Eqs.(2.4) and (2.5), depends on the suffix k [TT] . 
[20] . If p v in Eqs.(2.2) and (2.3) is replaced by p k , one would be led to 

(x u0 ,t u0 ) = Nl y^U^kUakU^Uai exp{i(p [H] x u0 - E [kl] t vQ )} 

k,l 

exp{--al(x u0 - v k t uQ ) 2 - -crl(x u0 - vit u0 ) 2 } 

= N l u ^U ak U lx iU a i cos(p[ k i]X u0 - EmUo) 
k,l 

exp{-al(x u0 - v {k i)tuo) 2 }(ixp{-^(rl(v [kl] ) 2 tl }, (A.l) 
which may be approximated as 



k,l 

1 

T 



cos{(p[ W ]i;(fci) - E [kl] )t vQ } exp{-^a^(% z] ) 2 ^ }, (A.2) 



where p\ k n = p k — p\. With a little algebra, one may verify that 

P[kl]V(ki) ~ E [kl] — — h [kl] . (A.S) 

In practical experiments, D = \x v q\ and t u o are not measured in an independent 
way; instead, the former is measured, while the latter is inferred from the former. 
Such a situation corresponds to t v or t to be integrated out in Eq.(A.l). One then 
obtains 



I/O; 



2E\ 



exp{-^(^Jo - ~'XTT )} exp{- }. (A.4) 
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Equations (A. 2) ~ (A. 4) coincide in the one-dimensional case with the corresponding 
equations derived in [TT] and [2U])@ It is interesting to observe the following. 

1. In Eq.(A.2), there appears a correction term to the "standard oscillation pe- 
riod", 4irE^i)/\ml — mf\, which appears in Eq.(2.10) or Eq.(2.13). Although 
this correction term vanishes in either of the equal-energy, equal-momentum 
and equal-velocity cases, it is rather artificial, in the framework of the wave- 
packet treatment presented in this paper, to assume Ek = Ei or v k = v\. 

2. In Eq.(A.4), each of oscillating factors is multiplied by two suppression factors: 
one which corresponds to that already present in Eq.(A.2), and in Eq.(2.10) or 
Eq.(2.13) as well, and the other which gives rise to another necessary condition 
for neutrino oscllation to be significant: \E[ k i]\/ a/2(£| + vf) < er„. 

Appendix B. Some algebra relevant to Sec. 4 

If Equation (4.1) is substituted, (v)k,i and (X)k t i (defined by Eq.(3.11)) and (u,w)k,i 
(given by Eq.(3.10)) may be expressed as 

(v)k,i = (v)(ki) ± 2^V[ k i], (X)k,i = (A) ( h) =F 2^V[ k i]t u , 
(u,w)k,i = (u, w)(ki) ± (u,w) kl + (u,w)h, (B.l) 

where 

(u, w) {M ) = — {^o-J«[(fcO^]^[(fcOM] + a l a l^m)n]W[(ki)7T] + crlalu^ n] w [fin] }, 

(«> w)'ki = ^ a l{ a l(u[M]™WM + tiw)fi]W[ki]) + (rl{u [kl] w mn] + u m7r] w [H] )}, 

w = ^^l{ a l + ^l)^[ki]W[kiy (B.2) 

Accordingly, a k j, bkj and c^i may be expressed as 

a-k,i = ct(ki) ± a k i + a 'L 
h,i = b( H ) ± b' M + b'l t , 

c k> i = c( H ) ± c kl + cl u (B.3) 



20 In the case in which p k and x u q may be expressed as pk — Pv\pk\l\pv\ (for any k) and as 
$ v q = p v \x vn \j\p v \, one has 

- rp ™ 2 -m 2 (jpfcj - \pi\)(\v k \ - \vi\) F 

P[kl]V(kl) ~ k[kl] = op JT, &lkl], 

P[kl] - -a , -Q &[kl] - \Vk\ - \Pl\ - ,-. | 2 . iQ -^[M] 
V% + Vf 1 J IVk] 2 + \Vi\ 2 



mj-mf (\vk\-\vimk\vi\~Ej\vk\) 
\Pk\ + \pi\ {\Pk\ + \pi\)(\vk\ 2 + \vi\ 2 ) 



E 



[kl] , 



2 2{v {kl) x u0 ) 2 {{\y k \ + I^DI^qI) 2 Jjvk\-\vi\f_ 2 

v^+v 2 M 2(\v k \ 2 + \v l \ 2 ) 2(\v k \ 2 + \v l \ 2 ) lXv01 ■ 
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where 

«(*!) = (v,v)(ki), a'ki = (v,v)'m, a'a = 

b m = (v,X) {kl) , b' kl = (v,X)' kl , b" kl = (€,jf)» M , 

c ikl) = (X,X) (kl) , c' kl = (X,X)' kl , d' kl = (X,X)l t . (B.4) 

Substituting Eq.(B.3) into Eq.(4.2), which defines ( k [, one may derive Eq.(4.6). On 
the other hand, substituting Eqs.(B.2) and (B.4) into q^) ~ a {ki)t\ k i)i one ma y derive 



c ( H) - a m t 2 {kl) = J Yl 



8a s a( k i\ 



(J K (T X (J K'^ \'i. V [K\]iX[ K '\']j - V[ K >x']jX [K xi 

Equation (4.5) therefore implies Eq.(4.7). From this, it follows that0 



b(ki) = —a(ki)to, C(fci) = d( k i)tl, 

Ki — —^kiitv + c 'ki — a kfivtoi 

t k>l = l(t v + t )- ^I'/i oll) (U-to). (B.6) 
2 2(a {H) ± a kl + a kl ) 

With the aid of Eqs.(B.5) and (B.6), one may verify that ( k i, Eq.(4.6), reduces to 
Ckiitu - W with C« given by Eq.(4.8)@ 
Similarly, if Equation (4.1) and 

Ek,i = E(ki) ± -Em (B.7) 

are substituted, AE k ,i and AE k j (defined respectively by Eq.(3.9) and Eq.(3.14)) 
may be expressed as 

1 ~ ~ , 1 ~, 

AE k j = AE( k i) ± -E[ H ], AE kj i = AE( kV) ± -AE W (B.8) 



where 

AE' kl = E [kl] - ^v [H] Ap. (B.9) 



(7 Z 



Substituting Eqs.(B.3) and (B.8) into Eq.(4.12), which defines rj kh one readily ver- 
ifies Eq.(4.16). Substituting Eqs.(3.9), (B.l) and (B.8) into Eq.(3.20), one may 
derive Eq.(4.20), which, if Equations (B.5) and (B.6) are further substituted, gives 
Eq.(4.21). 



Note also that the following relations hold independently of Eq.(4.7): 

ill If , If If ,2 m 

b ki = -a k iU, c kl = a kl t v . (B.5) 
It is not difficult to verify that Ckl > 0- 
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